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We study an elementary two-player card game where in each round players compare cards and 
the holder of the smallest card wins. Using the rate equations approach, we treat the stochastic 
version of the game in which cards are drawn randomly. We obtain an exact solution for arbitrary 
initial conditions. In general, the game approaches a steady state where the card densities of the 
two players are proportional to each other. The leading small size behavior of the initial card 
densities determines the corresponding proportionality constant, while the next correction governs 
the asymptotic time dependence. The relaxation towards the steady state exhibits a rich behavior, 
e.g., it may be algebraically slow or exponentially fast. Moreover, in ruin situations where one player 
eventually wins all cards, the game may even end in a finite time. 

PACS numbers: 02.50.Ey, 05.40.-a, 51.10.+y, 89.65.-s 



Numerous phenomena in social and economic sciences 
involve multiple interacting agents. The interaction be- 
tween these agents often leads to exchange of quanti- 
ties such as capital, goods, political opinions, etc. 
Games are widely employed in modelling collective be- 
havior especially in the context of economics jf|, with 
recent examples ranging from evolution of trading strate- 
gies in a stock market 0-^] to bidding in auctions |l0| . 
Games can often be regarded as many body exchange 
processes resembling collision processes [ fU| |, and there- 
fore their dynamics may be described by suitably adapted 
kinetic theories [Q. Here, we investigate a stochastic 
null strategy card game. By considering the "thermody- 
namic limit" where the initial number of cards is infinite, 
we show that rate equations provide a natural framework 
for analyzing game dynamics. 

Our two-player game is defined as follows. Each player 
starts with a certain number of cards. At each round 
players draw a card randomly from their deck and com- 
pare the card values. The player with the smaller card 
wins the round and gets both cards. This is repeated 
ad infinitum or until one of the players gains all cards. 
This stochastic adaptation of the elementary card game 
"war" is also motivated by a recently introduced auction 
bidding model where the lowest bid is rewarded JIo| . 

Our main result is that one specific aspect of the ini- 
tial card distribution, namely, the small size tail, gov- 
erns the dynamics of the game. Let us denote by A 
and B the two players, and let their initial card den- 
sities be ao(x) and bo(x), respectively. In the long time 
limit, a steady state is approached with the card den- 
sities of both players being equal to a fraction of the 
total card density, a^x) = a[ao(x) + bo(x)], with a oc (x) 
the limiting card density of player A. While a family 
of steady state solutions characterized by the parameter 
< a < 1 is in principle possible, the leading small size 
behavior of the initial distributions selects a specific value 
a = linu^n — r^rfu ? ^ ■ Moreover, the next leading cor- 

a (x)+b a (x) ' o 



rection determines how the system approaches the steady 
state. The corresponding time dependent behavior may 
be algebraic or exponential. Interesting behaviors also 
occur when one player captures all cards. In this case, 
the game duration maybe finite or infinite. Additionally, 
using numerical simulations wc show that the theoreti- 
cal predictions concerning the game duration extend to 
deterministic realizations of the game. 

Let the initial number of cards of player A and B be 
Na and Nb, respectively, and let the total number be 
N = Na + Nb ■ We shall take the thermodynamic limit 
Na,Nb,N — > oo such that the initial fractions Na/N 
and Nb /N are fixed. Then the card exchange process can 
be conveniently described via rate equations for a(x, t) 
and b(x, t), the densities of cards with value x at time t 
for players A and B, respectively. These densities evolve 
according to the nonlinear integro-differential equations 
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—a(x, t) = R(x, t), —b(x, t) = -R(x, t), (1) 
with the gain (loss) term R(x, t) given by 
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Here 
A{t) = 



dx a(x, t), 



B(t) = / 
Jo 



dxb(x,t) (2) 



are the fraction of cards possessed by players A and B, 
respectively. Clearly, 



A{t)+B(t) = 1. 



(3) 



The rate equations ([j]) reflect the nature of the game as 
the rate by which player A gains (loses) cards of value x is 
proportional to the fraction of his opponent's cards which 
is larger (smaller) than x. The overall factor L4B]" 1 en- 
sures that on average, every opposing pair of cards comes 
into play once per unit time. The minimal card value was 
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tacitly set to zero as the process is invariant under the 
transformation x — > x + const. 

Besides the obvious conservation law (^), two other 
conservation laws under ly the process. First, the total 
number of cards of a given value is conserved, 



a(x, t) + b(x, t) — uq(x), 



(4) 



where uq(x) — ao(x) + bo(x) is the initial total density. 
Second, the density of the minimal card remains con- 
stant throughout the evolution: a(0, t) — ao(0), where 
ao(x) = a(x,t = 0), and similarly for B. 

The steady state behavior is quite generic. The cumu- 
lative card densities, A(x, t) — J Q dya(y, t) and B(x, t) — 
Jq dy b(y, t) 7 satisfy A' /A = B'/B in the long time limit. 
Thus, Aoo(x) oc Boo(x), and consequently, the limiting 
card densities, a^x) = A'^x) and boo(x) = B' oc (x), 
are proportional to each other. The conservation law (Q) 
implies that each of the limiting card densities equals a 
fraction of the overall card density 



a-oo(x) = au (x), boo(x) = (1 — a)uo(x). 



(5) 



In principle, for a given total card density uo(x), there 
is a family of steady state solutions characterized by 
< a < 1. Moreover, initial conditions where the densi- 
ties are proportional to each other do not evolve further 
in time regardless of a. Still, for a given initial condition 
a specific value of a is selected. This value can easily be 
found for a class of initial conditions with non- vanishing 
minimal card densities, uo(0) > 0. Consider the den- 
sity of the smallest cards x — 0. Equation (|B|) gives 
floc(O) = auo(0), while the second conservation law im- 
plies aoo(0) = a (0), and hence a = ao(0)/[ao(0) + &o(0)]- 
This simple argument suggests that the smallest cards 
govern the outcome of the game. In the following, we 
solve for the full time dependent behavior and show that 
in general, the limiting small size tail of the two distri- 
butions dictates a. 

To solve the time dependent behavior, we make two 
simplifying transformations. First, the overall rate by 
which the exchange occurs L4.B] -1 can be absorbed into 
a modified time variable r, defined via 



ds [A{s)B(s)Y 
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The second transformation essentially reduces any total 
density uq(x) to a uniform density by introducing the 
variable £ 



dy u (y). 



(7) 



Suppressing the explicit time dependence, the trans- 
formed card densities are found from the relations 
d£ — a(x)dx and b(^)d^ = b(x)dx. Clearly, these 
satisfy a(£) = a(x)/uo(x) and &(£) = b(x)/uo(x). In the 



following, we shall omit the bar. The conservation law 
(0) becomes 



a(£,r)+&(£,r) = l, 



(8) 



i.e., the transformed total density is uniform on the in- 
terval [0,1] (note that Eqs. (|) and (|) imply < £ < 1). 

The above transformations simplify the evolution 
equations, and given the linear dependence (||), it suf- 
fices to solve for a 
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dV a{V, r) - a(£, r) / dr?6(r?,r). 
Jo 



Replacing &(£, r) with 1 — a(£, r) linearizes this equa- 
tion Jpa(£, r) = J^dr]a(ri,T) — £a(£,r), and differentiat- 



ing with respect to £ yields further simplification 



(9) 



Integrating over r and then over £ we arrive at our pri- 
mary result, the exact time dependent solution for arbi- 
trary initial conditions: 



d V a' ( V )e-^. 



(10) 
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Hereinafter we utilize the notations ao(£) = a(£,' 
a o(0 = l£ a o(0, and a = a o(£ = 0). 

Let us again consider the steady state. In the long time 
limit t — ► oo, the integral in ( |ic| ) vanishes and the den- 
sities become uniform a(£,r) — > a and b(^,r) —> I — a. 
Hence in terms of the original variable x, both densi- 
ties are proportional to uq(x) according to Eq. (|b|), with 
a = a o(£ = 0) = ao(x = 0)/uq(x = 0). Even when 
uq(x) vanishes or diverges as x — > 0, the parameter a is 
well-defined and using l'Hopital rule, it is given by 



,. a (x) 
urn — — ■ . 

x-,0 a (x) + b Q (x) 



(11) 



Thus, if the two distributions exhibit different leading 
behaviors, say lim^^o (%) / a o (%) — 0, then player A 
eventually ruins player B. Hence, the small card tail 
x — > provides the necessary selection criteria deter- 
mining which of the family of solutions (|^) is eventually 
selected by the dynamics. 

We now study the approach to the steady state. For 
example, the density A(t) = f Q d£ a(£, r) is given by 
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d£(l-£K(0e 



(12) 



While the steady state behavior is determined by the 
leading small argument behavior of ao(£), the relaxation 
towards the final state is governed by the correction to 
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the leading behavior. Let us consider the following small 
argument behavior 

a (£)~a + 7 £ 5 C^O, (13) 

with 6 > 0. Then, one has A(t) - a ~ -fT(8 + 1)t~ s . 
However, in terms of the actual time variable t, a richer 
variety of behaviors is exhibited. 

First, suppose that the system approaches an active 
steady state, i.e., < a < 1. Then from Eq. @ we 
obtain t — > a(l — a)r, and therefore 

i^-a-Cr*, i^oo (14) 

with C = "fT(8 + 1) [ck(1 — a)] 5 - Hence, the approach is 
generally algebraic. 

Next, suppose that one player, say A, eventually ruins 
their opponent, i.e., a = 1. Then dt/dr ~ -B(t) ~ t^* 5 
and consequently, t ~ r 1_<5 . Therefore, for 5 < 1 rep- 
resenting weak initial advantage of the eventual winner, 
the game duration is infinite: 

i-A(t) ~ *<i; (15) 

I „— const xt (5=1 

In the complementary situation of strong initial advan- 
tage for the eventual winner, S > 1, the game duration is 
finite: 

A(i/) = 1. (16) 

The terminal time can be determined from the integral 
t f = J °° drA{r) [1 - A(t)}. Using Eq. @ and recalling 
that a = 1 yields this time as an explicit function of the 
initial conditions 

tf=-[dt^a' (0 (17) 
Jo 5 

JO JO ?1 + ?2 

For example, the initial density a (£) = 1 — £ 2 yields 
t f = + j| In 2 » 0.87269. Additionally, the time de- 
pendent approach towards the final state is algebraic, 

l-A(t)~(t/-t)A, (18) 

sufficiently close to the terminal time t tf. As ex- 
pected, the density decreases linearly with time when the 
disparity between the two players becomes very large in 
the limit S — ► oo. 

Thus if the system approaches a trivial steady state 
with one player winning all cards, the temporal behavior 
can be algebraically slow or exponentially fast. Moreover, 
every positive power can be realized. Remarkably, if the 
initial disparity between the two players is sufficiently 
large, the game ends in a finite time. Interestingly, such 
disparity is expressed only in terms of the density of the 



smallest cards, with the larger cards practically irrelevant 
to the game outcome. 

Next, we analyze the time dependent evolution of the 
entire card density, not simply the overall number den- 
sity. Evaluating the leading behavior of the density (|l(]) 
in the long time limit, we find that the density exhibits 
a boundary layer structure 

a ^)-^{% s+1)T s f«£; d9) 

The scale £o ~ t _1 underlies the distribution. Cards 
larger than this scale have already relaxed to the limit- 
ing distribution, while cards smaller than this scale have 
yet to exchange hands and hence, are still dominated by 
the initial distribution. In other words, smaller cards are 
slower to equilibrate, consistent with the fact that the 
small size tail dominates the asymptotic behavior. 

We now briefly discuss the case where the number of 
card flavors is finite, or in other words, discrete card dis- 
tributions 

k 

a(x, t) = y~] a n (t)S(x - x n ), (20) 

n=l 
k 

b(x, t) = b n (t)6(x - x n ), 

71=1 

with x\ — and x n < x n +\. The discrete version of the 
rate equations can be written and solved directly using 
a series of transformation which mimics the ones used 
above. Instead, we shall insert the initial conditions (|20| ) 
in the general continuous case solution (|l0|). 

Denote again u n (t) = u n (0) = a„(0) + 6„(0) the to- 
tal card concentration. The variable = Y^m=i u m(®) 
plays the role of £ and the time variable r remains as in 
Eq. (|). The solution © reads 

«n(r) = qi(0) / a m (0) _ a m ~i(0) \ imT (2 , 

Un{0) «l(0) ^V«m(0) "m-l(0)y ' 1 ' 

Since all terms in the summation eventually vanish, the 
two players approach a limiting distribution which is pro- 
portional to the initial distribution a n (oo) = au n (0) with 
a = ai(0)/ui(0), in accordance with Eq. (|ll|). In gen- 
eral, the approach to the steady state is exponential. We 
first discuss the case < a < 1. Since Aoc — a, one 
has t — » a(l — a)r asymptotically. Hence, the relaxation 
towards the steady state is exponential 

A{t) - a - e - constx ' . (22) 

In the complementary case when one player wins all 
cards, a = 1, the approach is dominated by the first non- 
vanishing term in the summation, namely, the first non- 
vanishing b n (0). In this case, dt/dr cx exp(— const x r), 
and consequently, the game duration is finite as in 
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Eq. (|l(|). In summary, the behavior in the discrete case 
is different from the continuous case in that the time de- 
pendent behavior is generally exponential. An additional 
difference is that when one player captures all cards, the 
game duration is always finite. 

In the above, we discussed games with an infinite num- 
ber of cards. Nevertheless, one can apply the above re- 
sults to realistic situations when both players start with 
a finite number, say N, cards. We note that the time unit 
used above is of the order iV 2 rounds in an actual game. 
For the case 5 > 1 one therefore predicts a duration 

Tf ~ TV 2 , (23) 

with Tf the number of rounds. The duration in the 
marginal 5=1 case can be estimated using the aver- 
age time it takes for player B to get down to one card 
B(t) — N^ 1 . Utilizing the exponential decay of B(t), we 
find that there is an additional logarithmic dependence, 
T f ~ N 2 InJV, in this case. 
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FIG. 1. Duration of the game as a function of the number 
of cards. Shown is Tf the number of rounds a game lasts on 
average versus N, the initial number of cards. The results 
represent an average over 10 4 realizations. A line of slope 2 
is plotted for reference. 

Monte Carlo simulations are consistent with these pre- 
dictions. In the simulations, each player starts with TV 
cards drawn from a uniform distribution in the range 
< x < 1. Eventually, the player holding the smallest 
card wins. Our theory describes the stochastic realization 
of the game where cards are drawn randomly from the 
deck. We also examined the deterministic case where the 
card order is fixed throughout the game. In this version, 
the winner of a round places both cards in the bottom 
of the deck. In both cases, we find diffusive terminal 
times as in Eq. (^). Nevertheless, the two cases differ 



with the stochastic game ending faster than the deter- 
ministic one (see Fig. 1). Additionally, we find that this 
diffusive time scale characterizes the entire distribution 
of terminal times, as fluctuations in the terminal time are 
proportional to the mean {Tf) — {Tf) 2 oc {Tf) 2 . 

In closing, we studied a stochastic two-player card 
game using the rate equations approach. We found that 
extremal characteristics of the initial conditions select a 
particular steady state out of a family of possible solu- 
tions. Eventually, the card densities of the players be- 
come proportional to each other. However, the players 
generally possess different overall number of cards and it 
is even possible that one player gains all cards. The ap- 
proach towards the steady state exhibits rich behavior. 
Large cards tend to equilibrate faster than small cards, 
and the distribution develops a boundary layer structure. 
The time dependent behavior of the overall density is al- 
gebraic in cases where a steady state is approached. In 
the complementary case where one player gains all cards, 
the game may end in a finite or an infinite time. The rela- 
tive initial advantage of the winner, characterized by the 
correction to the leading extremal behavior, determines 
the game duration in this case. 
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